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ABSTRACT 

We analyze the ROSAT Deep Cluster Survey (RDCS) to derive cosmological constraints from the 
evolution of the cluster X-ray luminosity distribution. The sample contains 103 galaxy clusters out 
to z ~ 0.85 and flux-Hmit f/„„ = 3 x 10"" erg s'^cm'^ (RDCS-3) in the [0.5-2.0] keV energy 
band, with a high-rcdshift extension containing four clusters at 0.90 < z < 1.26 and brighter than 
Fiim = 1 X 10^*^ erg s^^cm^^ (RDCS-1). We assume cosmological models to be specified by the mat- 
ter density parameter f2m, the r.m.s. fluctuation amplitude at the 8/i~*Mpc scale cts, and the shape 
parameter for the CDM-like power spectrum P. Model predictions for the cluster mass function are 
converted into the X-ray luminosity function in two steps. First we convert mass into intra-cluster gas 
temperature by assuming hydrostatic equilibrium. Then temperature is converted into X-ray luminosity 
by using the most recent data on the Lx-Tx relation for nearby and distant clusters. These include the 
Chandra data for seven distant clusters at 0.57 < z < 1.27. From RDCS-3 we find = 0.35t°:i^ and 
CTg = 0.66lg'Q5 for a spatially fiat Universe with cosmological constant, with no significant constraint on 
r (errors correspond to la confidence levels for three fitting parameters). Even accounting for both the- 
oretical and observational uncertainties in the mass-X-ray luminosity conversion, an Einstein-de-Sitter 
model is always excluded at far more than the 3(t level. We also show that the number of X-ray bright 
clusters in RDCS-1 at z > 0.9 are expected from the evolution inferred at z < 0.9 data. Subject headings: 

Cosmology: theory - dark matter - galaxies: clusters: general - X-rays: galaxies 



1 Introduction 

According to the standard picture of hierarchical clus- 
tering, galaxy clusters arise from the gravitational collapse 
of exceptionally high peaks of the primordial density per- 
turbations. Therefore, they probe the high-density tail 
of the distribution of the cosmic density field, usually as- 
sumed to be Gaussian, and their number density is ex- 
ponentially sensitive to the cosmological scenario (e.g.. 
White, Efstathiou & Frenk 1983). The low-redshift cluster 
abundance has been used over the last decade to measure 
the amplitude of density perturbations on ~ 10/i~^Mpc 
scales (here h is the Hubble constant in units of 100 km s^^ 
Mpc""'^), while the redshift evolution of the cluster abun- 
dance refiects the growth rate of density perturbations, 
i.e., primarily depends on the matter density parameter 
fim (e.g., Oukbir & Blanchard 1992; Eke et al. 1998). Al- 



though very powerful in principle, this cosmological test 
faces two main problems in practical applications. First, 
theoretical predictions always provide the number density 
of clusters of a given mass, while the mass itself is never the 
directly observed quantity. Therefore, suitable assump- 
tions are required to relate an observational quantity to the 
actual cluster mass. Secondly, a cluster sample is needed 
which spans a large z-baseline, and is based on model- 
independent selection criteria, so that the search volume 
and the number density associated with each cluster are 
uniquely specified. In this respect, X-ray observations 
provide a very efficient method to identify distant clus- 
ters down to a given X-ray fiux limit, and hence within 
a known survey volume for each luminosity, Lx- For this 
reason, most of the studies which have used clusters as 
cosmological probes in the literature so far are based on 
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X-ray selected samples. Following the pioneering work 
based on the Einstein Extended Medium Sensitivity Sur- 
vey (EMSS, Gioia et al. 1990; Henry et al. 1992), deep 
imaging data from the ROSAT archive have been the basis 
for several serendipitous, flux-limited searches for high- 
redshift clusters (RDCS by Rosati et al. 1995, 1998; 160 
deg2 Survey by Vikhlinin et al. 1998; SHARC by Romer 
et al. 2000; WARPS by Jones et al. 1998; NEP by Gioia 
et al. 2001; see also Rosati et al. 2000 and Gioia 2000 for 
recent reviews). 

Estimates of the X-ray temperature, Tx, for subsets 
of these samples have opened the way to measure cluster 
masses to a fairly good, 15%, precision (e.g., Evrard, 
Metzler & Navarro 1996). The resulting X-ray tempera- 
ture functions (XTF) have been presented for both nearby 
(e.g., Henry & Arnaud 1991; Markevitch 1998; see Pier- 
paoli, Scott & White 2001, for a recent review) and distant 
clusters (e.g.. Eke et al. 1998; Donahue & Voit 1999; Henry 
2000), and have been compared with predictions from cos- 
mological models. The mild evolution of the XTF has 
been interpreted as a strong indication for a low-density 
Universe, with 0.2^ Jl,„^ 0.6. However, uncertainties re- 
lated to the limited amount of high z data and to the lack 
of a homogeneous sample selection for local and distant 
clusters could substantially weaken this conclusion (Co- 
lafrancesco, Mazzotta & Vittorio 1997; Viana & Liddle 
1999; Blanchard et al. 2000). 

An alternative method to estimate cluster masses is 
based on applying the virial theorem to internal velocity 
dispersions, cr,, , as traced by redshifts of member galaxies. 
This method leads to a rather precise determination of 
the mass function of nearby clusters (Girardi et al. 1998), 
although it is observationally very time consuming. The 
only statistically well-defined sample of distant clusters 
where this method has been applied is the CNOC survey 
by Carlberg et al. (1997), which contains clusters selected 
from the EMSS. Even in this case, the limited size and 
redshift extension of CNOC did not allow to stringent con- 
straints to be placed on from the distribution of cluster 
masses (Borgani et al. 1999b). 

It is worth noting also that the analyses realized so far 
on the evolution of the XTF and of the (Tt,-distribution 
combine samples of nearby and distant cluster, which have 
different selection criteria. In principle, this could compli- 
cate the comparison bc;tween low- and high-redshift data 
when establishing the evolution of the cluster population. 

A further method to trace the evolution of the clus- 
ter number density is to follow the evolution of the X- 
ray luminosity function (XLF). The relation between the 
observed Lx and the cluster virial mass is affected by 
the thcrmodynamical status of the intra cluster medium 
(ICM). Recent observations (Ponman, Cannon & Navarro 
1999) show an excess entropy in the ICM, not explained 
by gravitational processes. This demonstrates that non- 
gravitational heating and, possibly, radiative cooling sig- 
nificantly affects the Lx-M relation. Despite such com- 
plexities of the ICM physics, the X-ray luminosity has 
been shown to be a fairly robust diagnostic of cluster 



masses (e.g., Borgani & Guzzo 2001). Furthermore, the 
most recent flux-limited cluster samples contain now a 
large (~ 100) number of objects, which are homogeneously 
identified over a broad redshift baseline, out to z ~ 1.3. 
This provides a reliable way of combining data on nearby 
and distant clusters and a straightforward estimate of the 
selection function. 

Kitayama & Suto (1997) and Mathiesen & Evrard 
(1998) analyzed the number counts from different X-ray 
flux limited cluster surveys and found that resulting con- 
straints on flm are rather sensitive to the evolution of 
the mass luminosity relation. Sadat, Blanchard & Ouk- 
bir (1998) and Reichart et al. (1999) analyzed the EMSS 
and found results to be consistent with Q„i = 1- In our 
previous paper (Borgani et al. 1999, BRTN hereafter), 
we analyzed the XLF from the ROSAT Deep Cluster Sur- 
vey (RDCS), as derived by Rosati et al. (1998) at different 
redshift intervals, from a cluster sample shallower than the 
one analyzed here (sec Section 2.1, below). Wc wore able 
to set only moderate constraints on the density parameter; 
we found fl„i — 0.4 ±0.3 (90% confidence level) for a non 
evolving Lx T relation, with a critical density model still 
allowed by a moderate positive evolution of this relation. 
The weakness of these constraints was partly due to the 
method of analysis, and partly to the smaller size of the 
analyzed sample. More recently, a consistent result has 
been also found by Evrard et al. (2001), who compared 
the RDCS redshift distribution to results from the Hubble 
volume simulations. 

In this paper, we present the results on cosmological 
constraints from the analysis of the redshift and luminos- 
ity cluster distribution from the final version of the RDCS. 
This analysis differs from that in BRTN in several respects. 
First of all, the RDCS sample we analyze here is substan- 
tially enlarged (see Section 2.1), containing more than 100 
clusters selected to a fainter flux limit and out to a larger 
redshift, z< 1.3. Furthermore, in BRTN we fitted the XLF 
after computing it within finite Lx and z intervals. The 
analysis presented in this paper, instead, is based on a 
maximum-likelihood method which does not rely on any 
binning of the data and, therefore, has the advantage of 
exploiting all the information contained in the cluster dis- 
tribution within the whole portion of the {Lx,z) plane 
accessible to the RDCS. Finally, we base the luminosity- 
temperature conversion on the most recent observations, 
which probe now the whole redshift range sampled by 
RDCS clusters (e.g., Mushotzky & Scharf 1997; Donahue 
et al. 1999; Delia Ceca et al. 2000), with the notable ex- 
tension out to z ~ 1.3 from the Chandra observations of 
the ICM in the Lynx field (Stanford et al. 2001; Holden 
et al. 2001). We will devote particular care to verify- 
ing whether and by how much present uncertainties in the 
mass-luminosity relation weaken the derived constraints 
on the matter density parameter and on the amplitude of 
density perturbations at the cluster scale. 

The plan of the paper is as follows. After briefly intro- 
ducing the RDCS sample used for this analysis, we discuss 
the theoretical mass function and our approach to convert 
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observed X-ray fluxes into cluster masses. Finally, we de- 
scribe in detail the maximum-likelihood method that we 
apply to derive constraints on cosmological parameters. 
In Section 3 we present the results of our analysis and we 
discuss the main conclusions in Section 4. Unless other- 
wise stated, unabsorbed fluxes and luminosities refer to 
the [0.5-2.0] keV energy band. 

2 The RDCS analysis 
2.1 The sample 

The RDCS sample was constructed from a serendipitous 
search for extended sources in ROSAT PSPC pointings with 
exposure longer than 15 ksec. A wavelet-based algorithm 
was used to detect and measure the angular extent of X- 
ray sources. Over 160 cluster candidates were selected in 
180 PSPC fields as sources with an extent exceeding the 
local PSF with a 90% confidence level, which was statisti- 
cally determined by a control sample of several thousands 
sources (Rosati et al. 1995, 1998). 
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Fig. 1. — The cumulative redshift distribution, n(> z) of RDCS- 
3 clusters (continuous line) and of the high-z extension of the 
RDCS-1 clusters (dashed line). Also shown in the insert plot is 
the differential redshift distribution, N{z), for the same sam- 
ples. 

The RDCS sample that we consider in the following 
analysis is complete down to the flux limit = 3 x 

10"^'* ergs~^cm~^ (RDCS-3 hereafter) and contains 103 
spectroscopically confirmed clusters at z < 0.85 identi- 
fied over an area of approximately 50 deg^. This sample 
represents a substantial improvement with respect to that 
considered by BRTN, which included 70 clusters over an 
area of 32 deg^, above a flux-limit of 4x lO^^"* erg s~^cm~^. 
In Figure 1 we show the cumulative and differential red- 
shift distributions of the RDCS-3 clusters. RDCS-3 has 
overall a median redshift Zmcd = 0.29, with Zmax = 0.85; 
26 clusters lie at z > 0.5, and 4 clusters a,t z > 0.8. 
Several other clusters with 0.5 < z < 0.9 have been 
identified in the RDCS at Fx < 3 x lO^^"* crgs^^cm^^ 



these, however, do not belong to a complete sample (see 
below) and hence are not included in the present analy- 
sis. The largest flux is Fmax = 1-24 x 10^^^ crgs^^cm^^, 
with a median value of the flux distribution Fmcd = 
9.9 X 10"^"* ergs~^cm~^. Further details about the sam- 
ple, along with the presentation of the sky-coverage and 
an analysis of the XLF evolution are presented in a sep- 
arate paper (Rosati et al., in preparation). In addition, 
we also consider a deeper subsample of four clusters iden- 
tified down to Fiim = 1 X 10~^^ erg s~"'^cm~^ in the red- 
shift range 0.90 < z < 1.26 (RDCS-1 hereafter), whose 
z-distribution is also shown in Fig. 1. With the RDCS 
flux-limit and sky-coverage, we find, for a critical density 
Universe, that a cluster with Lx = 5 x 10''^'* ergs" "'^ has a 
searching volume Knax — 4.5 x 10'^( /i~^Mpc)^ at z > 0.5 
and Knax ^ 3.5 x W{h-^Mpc)^ for RDCS-3. 

As discussed in Rosati et al. (1998), the sky-coverage 
and sample completeness become uncertain at fluxes ^ 3 x 
10~^^ ergs~"'^cm~^. This is basically due to the fact 
that with ^ 50 counts, which roughly correspond to the 
above flux for the typical exposure time of the selected 
PSPC flelds, the detection and characterization of ex- 
tended sources becomes increasingly uncertain. Further- 
more, deep Chandra observations of the faintest RDCS 
clusters in the Lynx held (Stanford et al. 2001) show that, 
for such faint sources, PSPC-based A" -ray fluxes can be 
signiflcantly contaminated by emission from foreground or 
background point-like sources. For these reasons, we base 
our main analysis on the relatively bright sample, RDCS- 
3, while we will discuss whether the resulting constraints 
are consistent with the presence of the four clusters in the 
higher-z, lower-flux tail. Overall, our analysis will draw 
information on the evolution of the cluster number density 
from the widest redshift baseline presently accessible with 
deep cluster surveys. 

2.2 Modeling the Lx distribution 

Predictions for the Lx- and z-distributions of RDCS 
clusters are obtained by first computing the evolution of 
the cluster mass function and then by converting "theoret- 
ical" masses into observed X-ray fiuxes (see also Kitayama 
& Suto 1997; BRTN). 

The cluster mass function is usually written as 
n{M, z)dM = {p/M)f{v){dv/dM)dM, where p is the cos- 
mic mean density and v — Sc/(Jm{z). Here Sc is the critical 
density contrast for top-hat halo collapse, extrapolated at 
the present time by linear theory [Sc = 1.686 for Clm — 1), 
and (Ta/(z) the r.m.s. value of a top-hat density fluctua- 
tion at the mass-scale M at redshift z. Recently, Sheth & 
Tormen (1999) have proposed for f{i/) the expression 
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where a = 0.707, C = 0.3222 and q = 0.3, and showed 
that it provides a good flt to the halo mass function from 
N-body simulations (e.g., Governato et al. 1999; Jenkins 
et al. 2001). In the above equation the normalization is 
determined by the requirement J f{v)dv = 1. Eq.(P re- 
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duces to the standard Press-Schechter (1974) recipe for 
a = 1, C = 1/2 and q = 0. We assume each cosmologi- 
cal model to be specified by the matter density parame- 
ter, flm, and by the CDM-Hke power spectrum (Bardeen 
et al. 1986), with profile given by the shape-parameter 
r and normalization by the r.m.s. fluctuation amplitude 
within a sphere of 8/i~^Mpc comoving radius, as. Unless 
otherwise specified, in the following we assume flat spa- 
tial geometry provided by a cosmological constant term, 
fl\ = 1 — r2„i, consistent with recent small scale mea- 
surements of CMB anisotropics (de Bernardis ct al. 2000; 
Hanany et al. 2000). The value of f^A is known to have 
a relatively small effect on the evolution of the mass func- 
tion, and is far from being constrained by current cluster 
samples. 

Comparing a theoretical mass function with the ob- 
served Lx distribution of the RDCS requires a suit- 
able method to convert masses into X ray luminosi- 
ties in the appropriate energy band, [0.5-2.0] keV. As 
a first step we convert mass into temperature by as- 
suming virialization, hydrostatic equilibrium, and isother- 
mal gas distribution, according to the relation fc^T = 
1.38/3-^ Mf/^ [Qrn^yiriz)]^^^ {I + z) kcV (e.g.. Eke et al. 
1998). Here we take 76% of the gas to be hydrogen, M15 
is the cluster virial mass in units of 10^^h~^ Mq, (3 the ra- 
tio between the kinetic energy of dark matter and the gas 
thermal energy (/3 = 1 would be expected for a perfectly 
thermalized gas) and lS.mr{z) the ratio between the aver- 
age density within the virial radius and the mean cosmic 
density at redshift z (A^j^. = 187r^ ~ 178 for = !)• 
Although the assumptions on which the above relation 
is based have been recently questioned (e.g., Voit 2000), 
such a simple approach has been shown to reproduce fairly 
well the results from hydrodynamical cluster simulations 
(e.g., Bryan & Norman 1998, and references therein), with 
1~ /?~ 1.5. We assume for reference the value /3 = 1.15 
found by the Santa Barbara Cluster Comparison project 
(Frenk et al. 1999). We note that the M-Tx relation 
can be affected by the thermodynamics of the ICM. For 
instance, in their model for non-gravitational gas heating, 
Tozzi & Norman (2001) predict masses which can differ 
by about 20-30% with respect to those obtained from the 
above M-T relation. Possible deviations with respect to 
this relation are also indicated by observational data, when 
considering the mass within the internal cluster region at 
overdensity 5 > 500 (Finoguenov, Reiprich & Bohringer 
2001). Overall, we expect non-gravitational pre-heating 
to introduce fairly moderate changes in the M-Tx rela- 
tion used in our analysis. In order to account for them, in 
the following we will show the effect of changing (3 on the 
final model constraints. Finally, we assume 15% cluster- 
to-chister scatter in converting temperature into mass, as 
suggested by numerical simulations (Mctzler, Evrard & 
Navarro 1996). 

As for the relation between temperature and bolometric 



luminosity, we take the phenomenological expression 

(2) 

In this expression, is a dimensionless quantity and 

cIl^z) the luminosity-distance at redshift z for a given cos- 
mology, so that we explicitly factorize the redshift depen- 
dence induced by changing the spatial geometry of the 
cosmological background with respect to the Einstcin- 
de-Sitter (EdS) model. Several independent analyses of 
nearby clusters with Tx-~ 1 keV consistently show that 
Lg — 3 as rather stable results and a ~ 2.5-3 (e.g.. White, 
Jones & Forman 1997; Wu, Xue & Fang 1999, and ref- 
erences therein), with a rather small scatter, ^ 30%, es- 
pecially once cooling flow effects are taken into account 
(e.g., Markevitch 1998; Allen & Fabian 1998; Arnaud & 
Evrard 1999). For cooler groups, ^ 1 keV, the Lboi-Tx 
relation steepen, with a slope a ^ 5 (e.g., Helsdon & Pon- 
man 2000). As for the redshift evolution, Mushotzky & 
Scharf (1997) found that data out to z ~ 0.4 are consis- 
tent with no evolution for an EdS model (i.e., A ~ 0), a 
result which is consistent also with more recent data on 
cluster temperatures out to z ~ 0.8 (Donahue et al. 1999; 
Delia Ceca et al. 2000; Henry 2000). 

In Table 1 we provide an update of the results available 
on temperature determinations of distant {z > 0.57) clus- 
ters from recent Chandra observations. The two z > 1 
clusters (Stanford ct al. 2001) and the z = 0.57 cluster 
(Holden et al. 2001) have been observed in the Lynx field 
with 190 ksec ACIS-I exposure. The cluster CDFS-Cll 
has been detected as an extended source in the 1 Msec 
observation of the Chandra Deep Field South (Giacconi 
et al. 2001; Tozzi et al. 2001), its flux and tempera- 
ture having been estimated from about 350 counts (further 
details are presented in Giacconi et al. in preparation). 
The cluster 1WGAJ1226.9 at z = 0.89 has been serendip- 
itously discovered in the WARPS by Ebeling et al. (2001) 
and, independently, by Cagnoni et al. (2001) in a ROSAT 
PSPC blank field. Its temperature has been measured by 
Cagnoni et al. with a 10 ksec ACIS-S observation. As for 
the cluster MS1137, it has been observed with a 120 ksec 
ACIS-S pointing. We analyzed the corresponding Chandra 
archival data by applying the same procedure described by 
Stanford et al. (2001). The bolometric luminosity quoted 
in Table 1 is computed within an aperture radius of 60" . 
The resulting temperature, T = 5.7^Q g kcV, turns out to 
be consistent with, although somewhat more precise than, 
that determined by Donahue et al. (1999) from ASCA 
data. Results for the cluster MS 1054 have been obtained 
by Jeltema et al. (2001) from a 90 ksec exposiue with 
Chandra ACIS-S. Also in this case, the temperature is 
consistent with, although slightly smaller than that based 
on ASCA observations. 

The results for these distant clusters, along with other 
Lx-T data for both distant and nearby clusters, are shown 
in Figure 2. Overall, the results reported here constrain 
the evolution of the L^oi-Tx relation over the largest red- 
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Fig. 2. — The luminosity-temperature relation for nearby and distant clusters in two different cosmologies. Values of Lx assume 
here h = 0.5 for the Hubble parameter. The nearby clusters analyzed by Markevitch (1998) and by Arnaud & Evrard (1999) are 
indicated with small triangles and squares, respectively. Large circles are for the compilation of clusters at 0.5^ zi$ 0.8 reported 
by Delia Ceca et al. (2000). The large squares are for the compilation of distant (0.57 < z < 1.27) clusters recently observed with 
Chandra and reported in Table 1. The dashed lines indicate the Lboi-Tx relation of eq.(|^ with Le = 3 and a = 3, for A = (lower 
lines) and A — 1 computed at z = 1 (upper lines). 



shift interval probed to date. We stress that the high red- 
shift points always lie below or on top of the local relation, 
with the possible exception of the cluster at z = 1.27, 
whose X-ray spatial distribution suggests this might not 
be a relaxed system. In general, these data demonstrate 
that it is reasonable to assume A < 1, i.e., at most a mild 
positive evolution of the Ltoi-Tx relation. 
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Redshift 
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A03 






RXJ0848+4452" 


1.26 


n rtq + 0.83 

•^■^'-'-0.46 


5.58 


f 1.52 
-0.85 


5.8 


+2.4 
-1.3 


RXJ0848+4453" 


1.27 




+0.46 
^•^'-0.29 


1.6 


+ 1.2 
-0.6 


RXJ0848+4456'' 


0.57 


1 iq + 0.34 

^■^■-'-0.34 


1 6<)+°-^^ 
^•"^-0.51 


3.6 


+0.5 
-0.5 


CDFS-C11'= 


0.73 




0.22 


f0.09 
-0.09 


1.2 


+0.5 
-0.4 


1WGAJ1226. 9+3332'' 


0.89 




75.6 


+7.6 
-7.6 


10. 


-+4.0 
-■-3.0 


MS1137.5+6625 


0.78 




21.3 


+2.6 
-2.6 


5.7 


+0.8 
-0.6 


MS1054.4-0321= 


0.83 


34.5+^3:^ 


56.9 


+5.3 
-5.8 


10. 


1+1.0 
^-1.0 



Luminosity and temperature measurements for distant (z > 0.57) 
clusters from recent Chandra observations. Column 1: cluster ID 
(" Stanford et al. 2001; Holden et al. 2001; " Giacconi et al. 
2001, in preparation; Redshift from Ebeling et al. 2001, Lx and 
T from Cagnoni et al. 2001; " Jeltema et al. 2001). Column 2: 
spectroscopic redshift; Columns 3-4: bolometric luminosity, assum- 
ing h = 0.5, for an Einstein-de-Sitter (EdS) cosmology and for a 
flat low— density model with Qm = 0.3 (A03). Column 5: X-ray gas 
temperature and corresponding Icr uncertainties. 

Besides the relevance for the evolution of the mass- 
luminosity conversion, these results have profound impli- 



cations on the physics for the ICM. For instance, the model 
with constant entropy predicts —0.7 < A < 0.7 depending 
on the level of the entropy itself and its evolution with 
cosmic time (Tozzi & Norman 2001). Values of A in this 
range are significantly lower than the evolution expected in 
the self-similar case {A = 1.5). Therefore, both the shape 
and the non-evolution of the Lboi~Tx relation are well ex- 
plained in models with substantial non-gravitational pre- 
heating. 

In the following we will assume a — 3 and A = as ref- 
erence values, while we will also show the effect of changing 
both such parameters. Bolometric and K corrections to 
the [0.5-2.0] keV observed band are computed by using a 
Raymond-Smith (1977) model with Z = 0.3 for the mean 
ICM metallicity. The global scatter in converting X-ray 
luminosity into mass is estimated by adding in quadra- 
ture a 15% scatter in the M-Tx and a 30% scatter in the 
Lboi-Tx conversion. Its effect is then included in our like- 
lihood analysis by convolving the model luminosity func- 
tion with a Gaussian function having an r.m.s. scatter 
of 34%. Finally, for a given Lx, the flux is computed as 
F = Lx/[4^di(z)]. 

2.3 The analysis method 

We apply a maximum-likelihood (ML) approach, to 
compare the cluster distribution in the flux-redshift, 
(F, z), plane with predictions from a given cosmological 
model for a sample having the same flux-limit and sky- 
coverage as the RDCS. As a first step, we divide the {F, z) 
plane into elements of size dF dz and compute the model 
probability 

\{F,z)dFdz = n[M{F),z]^^^fsky{F)dFdz 

(3) 
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of observing an RDCS cluster with flux F at redshift z. 
Here dV{z) is the comoving volume element in the red- 
shift interval [z,z + dz] and fsky{F) is the RDCS flux- 
dependent sky-coverage (Rosati et al. 1998). If the bin 
size is small enough, such probabilities are always much 
smaller than unity, then the likelihood function C of the 
observed cluster flux and redshift is defined as the product 
of the Poisson probabilities of observing exactly one clus- 
ter in dF dz at each of the {Fi,Zi) positions occupied by 
the RDCS clusters, and of the probabilities of observing 
zero clusters elsewhere: 

£ = Y[[\iF,,z,)dzdFe-^'-''-'^^'''''^] J|e-^(^-^^)^^'^^. 

(4) 

Here the indices i and j run over the occupied and empty 
elements of the (F, z) plane, respectively. If we define the 
quantity S = —2hiC and drop all the terms which do not 
depend on the cosmological model, it is 



S = -2^1n[A(F,,z,)] 



dz I dFX{F,z) 



(5) 

(e.g., Marshall et al. 1983). In the above equation, Zmax 
represents the highest redshift at which the cluster identi- 
fication algorithm, on which RDCS is based, successfully 
detects extended sources and which, in principle, does not 
coincide with the highest-2: cluster identified in the sur- 
vey. The (1 -I- z)* surface brightness dimming is largely 
responsible for this high redshift cutoff. Using simula- 
tions as those shown in Fig.l of Rosati et al. (1999), 
we found Zmax = 1.5; at these redshifts the RDCS sam- 
ple becomes surface brightness limited for clusters with 
Lx — 10''''' ergs~^. We verified that final results are al- 
most left unchanged by taking Zmax = 2. Finally, we need 
to account for non-negligible errors in cluster fiuxes, which 
range from about 5% up to 35%, with a typical value of 
about 15%. To this purpose, we take clusters not to be de- 
fined as points on the (F, z) plane. Each cluster is spread 
along the F-direction using a Gaussian smoothing with 
r.m.s. amplitude equal to the flux error, ep- Therefore, 
instead of having zero or unity weight for empty and occu- 
pied cells in eq.(||), the i-th term contributing to the sum 
is assigned the weight 



E 



1 



2^4,, 



cxp 



dF, 



(6) 



where dF is its cell flux-width and the sum is over all the 
clusters having redshift between Zi — dz/2 and Zi + dz/2. 

3 Results 

We derive constraints on cosmological parameters by 
searching for the absolute minimum of 5* in the three- 
dimensional (r,crg,r2,„) parameter space, and compute 
confldence regions by allowing for standard increments 
AS". Cosmological parameters are varied within the follow- 
ing ranges: 0.02 < T < 0.4, 0.4 < as < 1.4, 0.1 < n„, < 1- 



In Figure 3 we show the constraints on the (rj,„, erg) plane 
for different values of F. As already mentioned, we as- 
sume for this reference analysis P — 1.15, Lq — 3, a = 3 
and A = 0. The most striking result is that the fairly large 
statistics and the wide z-baseline provided by RDCS allow 
stringent constraints to be placed, with Jl^ = 1 always ex- 
cluded at much more than the 3a level. The trend toward 
smaller 0^ for larger Fs is due to the fact that shallower 
spectra produce a steeper and more rapidly evolving mass 
function. Therefore, smaller fl„i values are required to 
compensate for the more rapid evolution. Although this 
effect is not large enough to qualitatively affect the re- 
sulting constraints, it has nevertheless some effect on the 
values of the best-fltting parameters. Overall, we obtain 
the following constraints on the cosmological parameters: 



nr. 



0.35 



-(-0.13 
-0.10 



0-8 = 0.66 



-1-0.06 
-0.05 



(7) 



(errors are la confidence levels for three interesting param- 
eters). No significant constraints are found for F, meaning 
that the sampled Lx range does not correspond to a mass 
range large enough to probe the shape of the power spec- 
trum. Indeed, RDCS spans about 2.6 decades in Lx- For 
the adopted M-Lx conversion, this corresponds to about 
1.3 decades in mass and, therefore, to about 0.4 decades in 
physical scales. For comparison, constraints on the power- 
spectrum shape from the galaxy distribution are typically 
derived by sampling the galaxy clustering over about two 
decades in physical scales. If we assume F = 0.25, as sug- 
gested by results from galaxy clustering (e.g., Dodelson 
& Gaztafiaga 2000; cf. also Efstathiou & Moody 2000), 
then the resulting constraints are flm = 0.28 ± 0.05 and 
as — 0.69 ± 0.06, where errors are now for two significant 
parameters. By assuming instead open geometry with van- 
ishing cosmological constant, this constraint changes into 
= 0.35±0.06 and as = 0.60±0.03. This shows the ten- 
dency of fiat models to favor slightly smaller values of il„i , 
as a consequence of the larger linear perturbation growth 
rate in the presence of a cosmological constant term. 

Also shown in Fig. 3 are the model predictions for the 
number of clusters expected in RDCS-1. The different 
curves in each panel show the loci of the rJm-crg plane cor- 
responding to a fixed number of RDCS clusters expected 
at z > 0.9 and F > 1 x 10~'^ erg s^'cm^^. Quite remark- 
ably, the four clusters detected in RDCS-1 can always be 
produced by cosmological models which lie inside the la 
contours defined by the RDCS-3. On the one hand, this 
result implies that the XLF evolution traced by brighter 
clusters at z^ 0.85 also extends at fainter fluxes out to 
the highest redshift reached by RDCS. On the other hand, 
it suggests that RDCS is not signiflcantly affected by in- 
completeness or unaccounted systematics also at fluxes 
10-"ergs-icm-2 < Fx < 3 x lO"'-* erg s-'cm-^. It 
should be said, however, that given the relatively small 
survey volume at these low fluxes, such systematics have 
in general a small impact on observed distribution func- 
tions. 



6 



CO 



1 

0.8 
0.6 
0.4 



,1111111111111 

- r=o.io - 


1 1 1 1 1 1 1 1 1 1 1 1 1 
r=0.25 - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

^K^, r=o.40 - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 rv i I i-i.] 


1 1 1 1 1 1 1 1 rf -1 ^ 1 1 1 1 1 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 



m 



m 



Q 



m 



Fig. 3. — Confidence regions on the Qm-crs plane for different choices of the power-spectrum shape-parameter F. Here a = 3, 
^ = and l3 = 1.15 have been assumed for the mass-luminosity conversion. Contours are la, 2a and 3a c.l. for two interesting 
parameters. The dotted curves show the the Qm-as relations corresponding to a given number of clusters, Nc, expected in RDCS-1 
at 2 > 0.9; Nc = 0.1, 1, 10, 30 from lower to upper curve. The solid curve correspond to the actually observed 4 clusters. 



In Fig. 4 we show the effect of changing in different ways 
the mass-luminosity conversion, after fixing F — 0.25. The 
main result from these tests is that < 0.6 at least at the 
3(7 c.l. for any variation of the mass-luminosity conversion 
within reahstic observational and theoretical uncertainties. 
Looking at the details of the effects, a positive evolution 
of the Lboi-Tx relation, which is only marginally allowed 
by data (e.g., Donahue et al. 1999; Delia Ceca et al. 2000; 
Fairley et al. 2000; see also Fig. 2), turns into smaller 
masses of more distant clusters for a fixed Lx- This in- 
creases the amount of evolution inferred for the mass func- 
tion and allows for a slightly larger 57^. A shallower slope 
of the local Lboi-Tx relation provides relatively smaller 
masses for less luminous clusters. Since low Lx objects 
are sampled at low redshift, this has the effect of decreas- 
ing the amplitude of the mass function at low z, so as to 
decrease its evolution, thus implying a lower fim- We have 
also verified that the observed steepening at the scale of 
galaxy groups (e.g., Helsdon & Ponman 2000) has negligi- 
ble effect on our results, only a few intermediate-z RDCS 
clusters being faint enough to be classified as groups. Fi- 
nally, a larger P implies a larger cluster mass for a fixed 
Tx, so that a larger erg is required to match the amplitude 
of the cluster XLF. In turn, a larger as slows down the 
evolution of the model mass function, thus allowing for a 
larger fim. 

These constraints on Vim are consistent with those found 
from some analyses of the XTF evolution. Eke et al. 

(1998) combined the temperature data for 25 local clus- 
ters by Henry & Arnaud (1991) with the sample of 10 
EMSS clusters at 0.3 < z < 0.4 by Henry (1997) and 
found — 0.40 ± 0.25 at the la c.l. Donahue & Voit 

(1999) used the low-z sample by Markevitch (1998) and 
enlarged the high-z sample by adding further 5 clusters at 
0.50 < z < 0.83; for flat geometry, they constrained the 
density parameter to lie in the range — 0.3±0.1 at the 
la c.l. Henry (2000) combined the low-z data by Henry 



& Arnaud (1991) with ASCA temperatures for 15 EMSS 
clusters at 0.3 < z < 0.6, and found = 0.44 ± 0.12 at 
the la c.l. for one fitting parameter, in the case of flat 
geometry. 

As already mentioned, such constraints from the XTF 
could be weakened by the lack of a homogeneous sam- 
ple of cluster temperatures selected at low and high red- 
shift. For instance, Viana & Liddle (1999) used the same 
data set as Eke et al. (1998) and showed that uncertain- 
ties both in fitting local data and in the theoretical mod- 
elling could significantly change final results. They found 
rim — 0.75 as a preferred value, with a critical-density 
model acceptable at < 90% c.l. Blanchard et al. (2000) 
derived a local XTF for clusters from the XBACs (Ebeling 
et al. 1996). The higher power-spectrum normalization 
they derived from this sample turned into a slower XTF 
evolution and, therefore, a higher As a result, they 
derived flm — 0.9 ± 0.3 at the la c.l. We stress here that 
the possible ambiguity connected with the choice of the 
low-z sample and the combination of different data sets 
at different redshift ranges is not an issue in our analysis, 
owing to the uniform selection provided by the RDCS over 
its entire redshift range. 

Besides the RDCS, the EMSS is the only other X-ray 
flux-limited sample that has been used to date to derive 
cosmological constraints. Sadat et al. (1988) used the 
EMSS redshift distribution, while Reichart et al. (1999) 
followed the evolution of the Lx distribution. After us- 
ing a slowly evolving Lx-T relation, they found flm to be 
consistent with unity and, therefore, significantly higher 
than the values preferred by our RDCS analysis. 

As for the amplitude of the power spectrum, our best 
fitting value of ag is somewhat smaller than those indi- 
cated by other analyses. For instance, fixing F = 0.25 and 
rim — 0.3, we find from the reference analysis shown in 
Fig. 3 that as = 0.67 ± 0.06 at the 3a c.l. for one inter- 
esting parameter, while Eke et al. (1998) find as — 0.8 
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Fig. 4. — Effect of clianging the mass-luminosity relation. Solid contours are from assuming F = 0.25, a = 3, ^4 = and l3 = 1.15. 
Each panel corresponds to changing one of the parameters defining the mass-luminosity relation. Contours have the same meaning 
as in Fig. 3. 



with about 20% uncertainty. We verified that this differ- 
ence can be partly explained by the fact that we used here 
the mass function by Sheth & Tormen (1999). At the ef- 
fective mass scale probed by RDCS this mass function is 
somewhat larger than that by Press & Schechter (1974), 
thus requiring a lower power-spectrum normalization to 
match data. We repeated our analysis with the standard 
Press-Schechter recipe and found as = 0.72 for the same 
choice of ftm and F. Furthermore, we verified that, for the 
same parameter choice, increasing /3 to 1.25 is enough to 
increase trg by a further 10%. 

Finally, we point out that the constraints obtained with 
the present analysis are more stringent than those we de- 
rived in BRTN. In that paper we used the XLF, binned in 
luminosity and redshift, for a smaller version of the RDCS, 
and found constraints on Q„i which were quite dependent 
on the assumed evolution of the Lx-T relation. Assuming 
non-evolution for this relation, we derived ~ 0.4 ± 0.3 
at the 90% c.l. for flat models. The much more stringent 
constraints derived in the present analysis are due to two 
main reasons: the larger number of clusters in the RDCS, 
extending to higher redshifts and fainter fluxes; and the 
new analysis which extracts more information than was 
previously possible with a simple grouping of RDCS clus- 
ters into luminosity and redshift intervals. Indeed, the 
likelihood function defined by eq.(H) conveys all the in- 
formation provided by the cluster distribution within the 
portion of the {F, z) plane accessible to RDCS. 

4 Conclusions 

In this paper we analyzed the evolution of the cluster 
number density, as traced by the ROSAT Deep Cluster Sur- 
vey (RDCS, Rosati et al. 1998, 2000) out to z ~ 1.3, to 
derive constraints on the r2„i and erg cosmological param- 
eters. Our analysis was aimed at understanding whether 
uncertainties in the relation between cluster mass, Af, and 
X-ray luminosity, Lx , prevent us from drawing firm con- 



clusions. In principle, a major source of uncertainty is 
related to the evolution of the M-Tx~Lx relation. The 
most recent data shown in Fig. 2 allow us now to trace 
the Lboi-Tx relation out to z ~ 1.3 (Stanford ct al. 2001), 
without evidence of significant evolution. This result also 
agrees with predictions from semi-analytical models of the 
ICM aimed at explaining the excess entropy observed in 
cluster cores (e.g., Tozzi & Norman 2001, and references 
therein). These findings are now significantly reducing the 
uncertainties on the evolution of the X-ray properties of 
cluster. 

As a main result we find that, within both theoretical 
and observational uncertainties in the M-Lx relation, the 
density parameter is always constrained to lie in the range 
0.1^ n,„^ 0.6 at the 3ct c.l. This demonstrates that X- 
ray luminosities as a function of redshift from deep X-ray 
flux-limited cluster surveys are indeed powerful tools to 
probe cosmological scenarios. 

Serendipitous cluster searches with Chandra and XMM- 
Newton archival data will lead to larger distant cluster 
samples within the next few years, and provide temper- 
ature information for a substantial number of objects. 
Flux-limited surveys will open the way to quantify the 
evolution of the cluster XLF with unprecedented accuracy, 
while the availability of many more cluster temperatures, 
even for incomplete samples, will provide a precise calibra- 
tion of the Lboi-Tx relation and its evolution. 

With this perspective, the main limitation to further 
tightening constraints on cosmological parameters will 
come from our theoretical understanding of what a cluster 
actually is. While we beneflted from the recent results on 
the thermodynamics of the ICM for the present analysis, 
another significant source of uncertainty stems from the 
dynamical aspects of clusters of galaxies. For instance, 
the limited validity of the assumptions on which the rela- 
tion between X-ray temperature and virial mass is based 
(e.g., spherical collapse and hydrostatic equilibrium, see 
Voit 2000), makes the connection with model predictions 
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somewhat uncertain. Furthermore, as observations are 
reaching the first epoch of cluster assembly, treating them 
as dynamically relaxed and virialized systems is undoubt- 
edly an oversimplification. In fact, hierarchical clustering 
scenario predicts that a fraction, between 0.3 and 0.6, of 
the z — 1 population of groups and clusters are observed 
less than 1 Gyr after the last major merger event and, 
therefore, are likely to be in a state of non-equilibrium. 

Although such uncertainties have been so far of minor 
importance with respect to the paucity of observational 
data, a breakthrough is however needed in the quality of 
the theoretical framework if high-redshift clusters are to 
take part in the high-precision era of observational cosmol- 
ogy. In this respect, hydrodynamical cluster simulations 
designed to include all the relevant ICM physics will play 
a fundamental role in the years to come. 
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